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Generalized Phase Synchronization in unidirectionally coupled chaotic oscillators 

Dae-Sic Lee, Won-Ho KyeQ Sunghwan Rim, Tae-Yoon Kwon, and Chil-Min Kin{j] 

National Creative Research Initiative Center for Controlling Optical Chaos, Pai-Chai University, Daejeon 302-735, Korea 

We investigate phase synchronization between two identical or detuned response oscillators cou- 
pled to a slightly different drive oscillator. Our result is that phase synchronization can occur 
between response oscillators when they are driven by correlated (but not identical) inputs from the 
drive oscillator. We call this phenomenon Generalized Phase Synchronization (GPS) and clarify its 
characteristics using Lyapunov exponents and phase difference plots. 
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Synchronization has been of much interest since Huy- 
gens' first description of it in two pendulum clocks on 
a wall 1]. The report that synchronization can be ob- 
served even in chaotic systems gave new rise to scien- 
tific attention to the phenomenon 0, . Over the past 
decade, synchronization in coupled chaotic oscillators has 
been intensely investigated for the understanding of its 
fundamental role in coupled nonlinear systems and the 
possibility of applications in various fields 0, 0, El, 0, • 
What characterizes synchronization is the convergence of 
the distance between the state variables of drive and re- 
sponse systems to zero due to weak interaction. Several 
different types of synchronization i.e., phase synchroniza- 
tion (PS) |1 ED, El, lag synchronization (LS) [H, com- 
plete synchronization (CS) 0, and generalized synchro- 
nization (GS) 13] have been observed in coupled chaotic 
systems. 

While CS, PS, and LS are observed in identical or 
slightly detuned systems (due to parameter mismatch) 
Q El EJ GS is observed in coupled oscillators with dif- 
ferent dynamics [l3|. When chaotic signals of a drive 
oscillator are fed into response oscillators, above a criti- 
cal coupling the response oscillators lose their exponential 
instability in the transverse direction and their state vari- 
ables converge to the same value. This convergence is the 
main character of GS. Since the attractors of response os- 
cillators converge to the same image in the GS regime, GS 
implies the emergence of a functional relation between 
drive and response oscillators such that xi = H(x2) |l3|, 
where xi and X2 are the state vectors of drive and re- 
sponse oscillators, respectively. 

In mutually coupled chaotic oscillators, there have 
been extensive investigations on the whole synchroniza- 
tion phenomena (PS, LS and CS) and various transi- 
tion scenarios clarified 0, E3- Most of the investiga- 
tions in unidirectionally coupled chaotic oscillators have 
been concentrated on GS transition and its applications. 
Also, some trials have been made to unveil the relation 
between PS and GS: Parlitz et al. E3 experimentally 
studied PS in unidirectionally coupled analog computer 
and Zheng et al. [l5[ theoretically demonstrated that 



GS can be weaker than PS depending on parameter mis- 
mat chs. Nevertheless, there remain unclarified questions 
concerning to PS in unidirectionally coupled chaotic os- 
cillators. Particular questions are how PS phenomena es- 
tablished in the drive-response and the response-response 
systems differ from each other and how they develop to 
GS. 

In this report, we investigate transition to PS in uni- 
directionally coupled chaotic oscillators driven by two 
different types of chaotic signal. We demonstrate that 
PS in the response-response system is induced by PS in 
the drive-response system when the driving signals are 
identical. However, we find that when correlated but not 
identical driving signals are used, PS is established in the 
response-response system but not in the drive-response 
system. We call this special PS phenomenon generalized 
phase synchronization (GPS) and discuss its relation with 
GS. 

How to define the phase for a chaotic system is an 
important issue and an active field of investigation in 
nonlinear dynamics. So far several methods, e.g. using 
phase space projection 19. fl 7| . H ilbert transformation^, 
and wavelet transformation [l^ etc., have been suggested 
and there have been done extensive investigation based 
on these methods. The method using phase space projec- 
tion is the most simple one to define the phase in Rossler 
oscillator and it enables us to use analytic treatment in 
analyzing the phase dynamics. We follow this approach. 
So the phase is defined by the simple geometric function: 
0i = arctan(^/a^), where i = 1 for a drive oscillator and 
i = 2, 3 for response oscillators. 

To demonstrate the conventional PS phenomenon es- 
tablished in the drive-response, we consider the unidirec- 
tionally coupled Rossler oscillators with slight parameter 
mismatch (see the caption of Fig. 1). The identical signal 
yi from drive oscillator is fed into two responses and the 
phase difference between drive and response oscillators 
can be written as follows 0, [HI 1 
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where <j> lk = 8\ - 6 k and R k \Jx k + y k and k = 2, 3. 
Here, Aco is the frequency mismatch between drive and 
response oscillators and £(£) is the fast fluctuating term 
which plays the role of effective noise. It is known 



that the above dynamics is governed by Type-I inter- 
mittency in the presence of noise and that PS is es- 
tablished when the channel width is deeper than the 
maximum of the effective noise [ll[. Accordingly, 
we can estimate the onset point of PS by consider- 
ing the fixed point condition: = 0. It leads to 
the equation: <j)\k = arcsin(2Ao; J R/ c /ei?i) where we ig- 
nore the effective noise term £ because it is negligible 
in the PS regime |9|- Accordingly the onset solution of 
the fixed point is cj)* k = tt/2 when 2AcuRk/cRi = 1. 
Thus the critical coupling can be estimated by e ~ 2Aco 
since Ri/Rk ~ 1 between slightly detuned systems 
Then the critical coupling for PS can be estimated by 
e c = 2Au = 2(cv d - uo r ) 0.03. 

Figure 1 shows the stroboscopic phase trajectories and 
probability distributions P(Qi, 2 ) of chaotic oscillators 1 
and 2 just above the critical point with the reference 
oscillator (see Ref. [l8| for our definition). One can see 
that PS occurs between oscillators 1 and 3 (Fig. 1 (a) and 
(b)) as well as between oscillators 2 and 3 (Fig. 1 (c) and 
(d)) as oscillators 1 and 2 take the preference directions 
(i.e., (exp(z0)) = / exp(i0)P(0)d0 + 0). This implies 
that the rotational symmetry on the projective attractor 
is broken due to PS transition, which is the indisputable 
evidence of PS 

0, E3, EH Accordingly, we understand 
that PS established in the response-response system is 
due to PS in the drive-response system, i.e., above the 
critical coupling (\)\ 2 = and (\)\ 3 = imply 023 = in 
Eq. (1). In other words, PS in the drive-response system 
coincides with PS in the response-response system when 
identical driving signals are used. 

Next, we consider two response Rossler oscillators 
driven by the correlated signals, x\ and 2/1, instead of 
identical ones: 

yi = uj d xi + 0.1652/1, 

11 = 0.2 -Mi On - 10), (2) 

x 2 -u r y2 - z 2 + e(xi - x 2 ), 
y 2 = cv r x 2 + 0.1652/2, 

1 2 = 0.2 + z 2 (x 2 - 10), (3) 

x 3 -ov r y 3 - 23, 

ys = oj r x 3 + 0.1652/3 + e(yi - 2/3), 

z 3 = 0.2 + z 3 (x3-10), (4) 

where the correlated signal x\ and 2/1 of oscillator 1 are 
fed into oscillators 2 and 3, respectively. In real systems, 
noise and delay in propagating channel are unavoidable. 
Thus the above system models a real situation in which 
two response systems are driven by correlated signals, x 
and x' where x' is a distorted version of x. We propose 
the above system for studying PS in unidirectionally cou- 
pled systems and its relation to GS. 

The difference dynamics between two response oscilla- 
tors is given by AX = AAX + H(t) where AX = x 2 — x 3 , 
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FIG. 1: PS at e = 0.04 in coupled Rossler oscillators: x\ — 
-Wdyi - zi, yi = oudXi + 0.152/1, zi = 0.2 + zi(xi - 10), 

X2,3 = ~U r y2,3 - £2,3, 2/2,3 = U r X 2 ,3 + 0.1652/2,3 + e(?/i - 2/2,3), 

Z2,3 = 0.2 + ^2,3(^2,3 — 10), where — 1.015 and u r — 1.0. 
(a) stroboscopic phase trajectory (black dots) of oscillator 1 
with reference oscillator 3 [18|. Gray dots show the whole 
attractor of oscillator 1 without stroboscopic sampling, (c) 
stroboscopic phase trajectory (black dots) of oscillator 2 with 
reference oscillator 3. (b) and (d) probability distributions of 
(a) and (c), respectively. 
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FIG. 2: Two largest transverse Lyapunov exponents when 
ou d = 0.7 and uo r = 1.0. There are two transition points, A 
and B. 



A = ((0,-^,-1), {uj r , 0.165,0), (0,0,-10)), and S = 
diag(e(xi -x 2 ), - e(yi - 2/3), z 2 x 2 - z 3 x 3 ). By iterating 
this dynamics, we can find transverse Lyapunov expo- 
nents describing the relative motion of oscillators 2 and 
3. In Fig. 2, we see two transition points, A and B, 
which, as we will see below, correspond to two different 
types of PS: namely, GPS at A and the conventional PS 
at B. 

The phase difference between drive and response oscil- 
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FIG. 3: GPS in Eq. (2)-(4) at e = 0.2: (a) stroboscopic 
phase trajectory (black dots) of oscillator 1 with oscillator 3 
|18|: (c) oscillator 1 with oscillator 2; (e) oscillator 3 with 
oscillator 2; (b) (d) (f) probability distributions of (a), (c), 
and (e), respectively. 
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FIG. 4: Temporal behaviors of drive and response oscillators 
in the GPS regime when e = 0.2: (a) x\ and X2 (b) x\ and 
X3 (c) X2 and £3. 



lators is given by 

<j) lk = Acj - £(0i, 6 k ) sin (j) lk +77fc(0i,0fc), (5) 
where, Acj = — c<; r , 

B(0iA) = 5^-0.15008^! + ^), 

m(0i,e k ) = £ ( 2fc - 5 ) gi sin( g 1 + g fc) 

— (0.015 — e) sin 0& cos 0/~ 

+ (— -sin 0i - — sm0fc). 
itl it/c 

We can see a /c-dependent term in ^(0i, Ok) which is due 
to the driving by correlated signals from the drive oscil- 
lator in Eq. (2)- (4). We can obtain the critical point for 
PS transition which is e c = — 2Acj = 0.6, in accordance 
with the argument of the former case (below Eq. (1)). 
The drive and response oscillators develop to a PS state 
at this critical value and PS between oscillators 2 and 3 
is induced above this critical coupling. The critical value 
of the critical coupling agrees with that of the transition 
point indicated by point B (e ~ 0.55). Thus we under- 
stand that B corresponds to conventional PS transition 
point at which the three oscillators develop to PS simul- 
taneously. 

We need to inspect the phenomenon at reference point 
A, which is described as crossing to the negative value 
in one of the Lyapunov exponents. Fig. 3 shows the 



trajectories in phase space and probability distributions 
near reference point A. A PS state appears only in the 
response-response system (Fig. 3 (e) and (f)) without 
PS in the drive-response system (because probability dis- 
tributions are not localized in Fig. 3 (a)-(d)), which is 
different from the case of the reference point B. We call 
this phenomenon GPS on the analogy of GS in which two 
state variables, X2 and X3 coverge to the same value. In 
GPS the phases 02 and 03 are bounded by a constant. 
Figure 4 shows the temporal behaviors of each oscillator 
at the same coupling strength as that of Fig. 3. We note 
that PS is established in oscillator 2 and 3 as the phases 
are mostly matched in Fig. 4 (c), while phase slippings 
appear in oscillators 1 and 2 (Fig. 4 (a)) or oscillators 
1 and 3 (Fig. 4 (b)), intermittently. Finally, we remark 
that GPS is a novel phenomenon characterized by PS 
in the response-response system, and thus it is different 
from conventional PS|9|- 

In Fig. 5, we also observe that GPS appears when 
the responses are slightly detuned with frequencies of 
1.0±0.005, and it destabilizes for frequencies of l.OzbO.Ol. 
This implies that GPS is a real phenomenon that should 
be experimentally observable pjj . And the attractor de- 
formations observed in Fig. 3 (e) and Fig. 5 (e) seem to 
be originated from the driving signal of different natural 
frequency. The similar phenomenon is often observed in 
unidirectionally coupled systems with different dynamics 
[l3| . It was shown recently that the phase defined 

by geometrical function like ours and the phase based on 
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FIG. 5: Appearance of GPS when two response oscillators 
are slightly detuned. The natural frequencies of response os- 
cillators are w r = 1.005 (of oscillator 2) and w r — 0.995 (of 
oscillator 3), respectively. Others parameters are those of Fig. 
3. 



the Hilbert transformation practically coincide. So we 
think our result is independent of the method of defining 
the phase. 

In conclusion, we have studied PS in unidirectionally 
coupled chaotic systems with parameter mismatch. And 
we have focusedly clarified the relationship of PS phe- 
nomena in the drive-response and PS in the response- 
response systems. When the driving signals are identical, 
PS in the drive-response system corresponds to PS in the 
response-response system and the system develops to GS 
as the coupling strength increases: PS^GS. When the 
driving signals are correlated (but not identical), PS is 
established in the response-response system but not the 
drive-response system. We call this phenomenon GPS. 
The GPS state transits to PS when coupling strenght in- 
creases. The results are confirmed by the analysis of Lya- 
punov exponents, phase trajectories, and time series. We 
expect that the GPS concept could be used for analyz- 
ing weak interdependences of data coming from weakly 
correlated systems such as neuronal systems 0, cardiac 
oscillators [2(|, and ecological systems [2l| etc. 
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